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Abstract. Radio loud active galactic nuclei present a remarkable variety of signs indicating the presence of periodical processes 
possibly originating in binary systems of supermassive black holes, in which orbital motion and precession are ultimately 
responsible for the observed broad-band emission variations, as well as for the morphological and kinematic properties of the 
radio emission on parsec scales. This scenario, applied to the quasar 3C 345, explains the observed variations of radio and 
optical emission from the quasar, and reproduces the structural variations observed in the parsec-scale jet of this object. The 
binary system in 3C345 is described by two equal-mass black holes with masses of » 7.1 x IO'^Mq separated by a 0.33 pc 
and orbiting with a period ~ 480 yr The orbital motion induces a precession of the accretion disk around the primary black 
hole, with a period of « 2570 yr The jet plasma is described by a magnetized, relativistic electron-positron beam propagating 
inside a wider and slower electron-proton jet. The combination of Alfven wave perturbations of the beam, the orbital motion 
of the binary system and the precession of the accretion disk reproduces the variability of the optical flux and evolution of the 
radio structure in 3C 345. The timescale of quasi-periodic flaring activity in 3C 345 is consistent with typical disk instability 
timescales. The present model cannot rule out a small-mass orbiter crossing the accretion disk and causing quasi-periodic flares. 
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1. Introduction 

Most of the radio loud active galactic nuclei (AGN) exhibit 
emission and structural variability over the entire electromag- 
netic spectrum, on timescales ranging from several hours to 
several years and on linear scales of from several astronomical 
units to several hundreds of parsec s. In AGN with prominent 
relativistic jets, this variability appears to be related to the ob- 
served morphology and kinematics of the jet plasma on parsec 
scales (e.g., Zensus et al. 2002 . Jorstad et al. 2001.) . In a number 
of radio-loud AGN, enhanced emission regions (components) 
embedded in the jet move at superluminal speeds along heli- 
cal trajectories (Zensus, Krichbaum & Lobanov .199 5 1996 
Zensus [r997> . and the ejection of such components is often as- 
sociated with optical and/or y-ray outbursts. 

In a growing number of cases, explanation of the observed 
nuclear variability and structural changes in parsec-scale jets 
is linked to the presence of supermassive binary black hole 
(BBH) systems in the centers of radio loud AGN. The BBH 
model was originally formulated by Begelman, Blandford & 
Rees (1980^, and it was studied subsequently in a number 
of works (e.g., Polnarev & Rees 1994 Makino & Ebisuzaki 
[19961 Makino HWTl Ivanov, Papaloizou & Polnarev [D^ . 
The BBH model postulates that an accretion disk (AD) exists 
around at least one of the two black holes (typically, around 
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the more massive one), and the resulting variability and struc- 
tural changes are determined by the dynamic properties of the 
disk itself and the BBH-AD system (this includes the disk 
and black hole precession, orbital motion, passages of the sec- 
ondary component through the AD). 

The most celebrated examples of successful application of 
the BBH model include OJ 287 (Sillanpaa, Haarala & Valtonen 
1988, Valtaoja et al. 2000 1, 3C273 and M 87 (Kaastra & Roos 
D.992I, 1928+738 (R oos, K aastra & Hummel [T^ . MrkSOl 
(Rieger & Mannheiml2000l and PKS 0420-014 (Britzen et al. 
2001 1. The epochs and structure of major outbursts in OJ287 
have been suggested to result from passages of the secondary 
black hole of a binary pair in OJ 287 through the accretion disk 
around the primary during the orbit (Lehto & Valtonen 1 1996> . 
The component ejection epochs are correlated with the optical 
outbursts, and the component trajectories show evidence for a 
heUcal morphology of the jet in OJ 287 (Vicente, Chariot & Sol 
11996'(. Several jet components observed in 3C 345 also move 
along distinct helical paths (Steffen et al. 119951 Lobanov 1996, 
hereafter lL96l . Lobanov & Zensus (1999, hereafter [CT^ 
have shown that shocks are not likely to play a significant 
role in the dynamics and emission of parsec-scale regions in 
3C345. Recently, the presence of binary black hole systems 
(albeit with extremely short precession periods and small or- 
bital separations) has been suggested for 3C 279 (Abraham & 
Carrara[T998j, 3C 273 (Abraham & Romero [T999t and 3C 345 
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(Caproni & Abraham I2004> . from the observed periodic varia- 
tions of the component ejection angle. 

The observed helical trajectories of jet components are a 
strong indication of precession or perturbation of the jet flow. 
The dynamics and emission of such perturbed outflows have 
been explained in the framework of the two-fluid model (Sol, 
Pelletier & Asseo .1989. Hanasz & Sol .1996 J describing the 
structure and emission of the jet in terms of an ultra-relativistic 
electron-positron (e'^) beam with Fb ~ 10 surrounded by a 
slower, electron-proton (e^p) jet moving at a speed ySj ~ 0.4 c. 

Observational evidence for the two-fluid model is reviewed 
in Roland & Hetem (1996 1. The presence of two fluids in 
extragalactic outflows was first inferred in the large-scale jet 
of Cygnus A (Roland, Pefletier & Muxlow 1988 1. The two- 
fluid model is supported by y-ray observations of MeV sources 
(Roland & Hermsen [T995l Skibo, Dermer & Schlickeiserl^ i 
and polarized structures in the compact radio jet of 1055+018 
(Attridge, Roberts & Wardle[T999}. The X-ray and y-ray ob- 
served in Cen A may be formed by the beams (Marcowith, 
Henry & Pefletier 1995; Marcowith, Henri & Renaud 1998 1. 
Direct detections of fast and slow jet speeds have been reported 
for Cen A (Tingay et al. 1998^ and M 87 (Biretta, Sparks & 
Macchetto 1999 1. Subluminal speeds, possibly related to the 
e^p plasma, have been observed in the double-sided jets in 
3C338 (Feretti et al. 1993 1. The two-fluid scenario has been 
used to explain the morphology and velocity field in the large- 
scale jet in 3C 31 (Laing & Bridle l2002ll2004t . 

Orbital motion and disk precession in the BBH-AD system 
perturb (but not disrupt) the beam, which leads to the com- 
plex variability and kinematic patterns observed in the compact 
jets. On scales of up to several hundreds of parsecs, the e- beam 
is not disrupted by Langmuir turbulence (Sol, Pelletier & Asseo 
[T989t and Alfven and whistler waves (Pelletier & Sol [T992l 
Achatz & Schlickeiser' 1 993V No strong Kelvin-Helmholtz in- 
stability should form on these scales, despite the significant 
transverse stratification of the outflow (Hanasz & Sol 119961 . 
Kinematic and emission properties of perturbed beams have 
been calculated (Roland, Teyssier & Roos 1994.. Desringe & 
Fraix-Burnet 1 1 9971 . and it has been shown that the observed 
variability and non-linear trajectories of jet components may 
indeed be caused by a perturbed plasma. In a more gen- 
eral approach, solutions for the restricted three-body problem 
(Laskar 1990 Laskar & Robutel 7995^ have been apphed for 
deriving the parameters of the perturbed beam and precessing 
accretion disk in the BBH system in the quasar PKS 0420-014 
(Britzen et al. l200n . 

In this paper, the dynamics of the BBH-AD system and the 
two-fluid model are combined together in order to construct a 
single framework for explaining optical flaring activity, kine- 
matic properties and internal structure of the compact parsec- 
scale jets. The combined model is developed to describe the 
quasar 3C345 . It is presented in Sect. |2] Section |3l provides 
an overview of the basic properties of the optical and radio 
emission from 3C345. In sections |4j|5l the model is applied 
to describe the evolution of the optical emission and radio jet 
in 3C345 after a strong flare in 1992 that resulted in the ap- 
pearance of a new superluminal jet component CI. 



Throughout the paper, the Hubble constant Ho - 70 h km 
s^'Mpc"', deceleration parameter qo - 0.5, and negative 
definition of the spectral index, S oc v " are used. For 
3C345 (z = 0.595, Hewitt & Burbidge 1993 1, the adopted 
cosmological parameters correspond to a luminosity dis- 
tance Z)l - 1.99/z^'Gpc. The corresponding linear scale is 
3.79/1"' pcmas"'. A proper motion of 1 mas/year translates into 
an apparent speed of 19.7/i"'c. Subscripts "b", "j", "c", and 
"m" are introduced to identify quantities related to the beam, 
jet, moving plasma condensation, and magnetic field, respec- 
tively. 

2. The binary black hole model 

The binary black hole model provides a general formalism for 
calculating the emission and kinematics of an outflow originat- 
ing from a BBH-AD system. The two-fluid description of the 
outflow is adopted (Fig.[0, with the following assumptions: 

1 . The outflow consists of an e^p plasma (hereafter the jet) 
moving at mildly relativistic speed /3j < 0.4 c and an e"^ 
plasma (hereafter the beam) moving at highly relativistic 
speed (with corresponding Lorentz factors Fb ~ 10). Both 
beam and jet may be transversely stratified. 

2. The magnetic field lines are parallel to the flow in the beam 
and the mixing layer, and the field is toroidal in the jet 
(Fig ID- The energy densities of the beam, jet and magnetic 
field are related as follows: ej » Sb ~ fim- The magnetic 
field lines are perturbed, and the perturbation propagates 
downstream at the Alfven speed, Va. 

In the two-fluid model, the e p jet carries most of the mass 
and kinetic power released on kiloparsec scales (large scale 
jets, extended lobes, hot spots). The e"^ beam propagates in 
a channel inside the jet, and it is mainly responsible for the 
formation of superluminal jet components observed on parsec 
scales. The beam interior is assumed to be homogeneous, en- 
suring that the beam rotation does not affect the observed evo- 
lution. Assuming a pressure equilibrium between the beam and 
the jet, a typical size of the beam is s O.lRj, which is of the 
order of 0.1 pc (Pelletier & Roland 1989 1. For AGN at distances 
of the order of 1 Gpc, this translates into an angular dimension 
of ~30 microarcseconds, and the beam remains transversely 
unresolved for very long baseline interferometry (VLBI) ob- 
servations at centimeter wavelengths. 

2.1. Geometry of the model 

Basic geometry of the model is illustrated in Fig. |2l Two black 
holes Ml and M2 orbit each other in the plane OfO^). The center 
of the mass of the binary system is in O. An accretion disk 
around the black hole Mi is inclined at an angle Qp with respect 
to the orbital plane. The angle f2p is the opening half-angle 
of the precession cone. The disk precesses around the ^ axis, 
which is parallel to the ^-axis. The momentary direction of the 
beam axis is then given by the vector MiJ. The line of sight 
makes an angle t with the precession axis. 
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Fig. 1. Two-fluid model for relativistic outflows. A fast, rela- 
tivistic beam is surrounded by a slower, possibly thermal, 
e^p outflow, with a mixing layer forming between the beam 
and the jet. The magnetic field, B, is parallel to the flow in the 
beam and in the mixing layer, and the field is toroidal in the jet. 
The magnetic lines are perturbed and the perturbation propa- 
gates outwards at the Alfven speed, Va. 



For a typical parsec-scale jet, this implies «b ^ lOe" cm . The 
beam moves along magnetic field lines in the e jet if the 
magnetic field is larger than the critical value 

B„it = i^nn^i^m.c^yi^ (2) 

where ^mix ~ nj is the density of the plasma in the beam-jet 
mixing layer (Pelletier, Sol & Asseo 1988 1. The magnetic field 
in the jet rapidly becomes toroidal (Pelletier & Roland[l990j. 
The beam is not disrupted by Kelvin-Helmholtz perturbations 
if the magnetic field is larger than 

Bkh = (27rnbr^mec2)i/2 (3) 

(Pelletier & Roland 1989). A typical density ratio in parsec- 
scale jets is = (nb/«j) - 0.01, which implies Bcrit ~ ^kh ~ 
50 mG. Sol, Pelletier & Asseo ( I1989> have demonstrated that 
the beam does not suffer from inverse Compton and syn- 
chrotron losses on parsec scales. Relativistic bremsstrahlung, 
ionization and annihilation losses are also insignificant on these 
scales. Diffusion by turbulent Alfven waves is not expected for 
beams with Fb ~ 10. 




/ Observer 



Fig. 2. Geometry of the BBH model. Mi and M2 denote the 
locations of two black holes orbiting around each other in the 
plane C^^O^). The coordinate system is referenced to the center 
of mass of the binary system (O). The accretion disk around 
the black hole Mi inclined at an angle Qp with respect to the 
orbital plane, such that the | axis (parallel to the ^ axis) is the 
precession axis. The momentary direction of the beam axis is 
given by the vector MiJ. 



2.2. Stability and dissipation of tiie beam 

The beam exists inside the jet for as long as its density 

Hb >Sj/(r^mec2). (1) 



2.3. Kinematics of tlie perturbed beam 

The relativistic plasma injected into the beam in the vicinity 
of the black hole Mi moves at a speed characterized by the 
bulk Lorentz factor Fc. The condensation follows the magnetic 
field lines. The perturbation of the magnetic field propagates 
at the Alfven speed Va = B(47rmenb)"*'^. For highly energetic 
beams, the classical density po - m^n^, is substituted with the 
enthalpy yU = p + Pjc^, where P is the pressure, p = po(l +s/c^) 
is the relativistic density and s is the specific internal energy. 
Taking into account relativistic corrections, and recalling that 
the enthalpy becomes Fb(p + P/c^) in a plasma moving at a 
bulk speed ySb = (1 - 1 /F^)"-'', the resulting phase speed of the 
Alfven perturbation with respect to the beam plasma is VA.rei = 

1/(1 -t-Fb/y^). Thus, for a typical e"^ beam propagating at 
Fb S 10 inside a jet with B S B^nt, the relativistic Alfven speed 

is a: 0.1 C. 

The observed trajectory of the beam is affected primarily 
by the orbital motion in the BBH system and the precession 
of the AD around the primary BH. The jet structure on parsec 
scales is mainly determined by the precession. The orbital mo- 
tion introduces a small 0.1 mas) oscillation of the footpoint 
of the jet, and it is not immediately detectable in VLBI images 
with a typical resolution of » 1 mas. Consider first the effect 
of the precession of the accretion disk. The coordinates of the 
component moving in the perturbed beam are given by 



iX = R(0 cos(d^pf - fep^ + (/>p) 
I = R(^) sin(wpf - kp^ + 4>p) 
,1 = ^ = ^(0, 



(4) 



where R(^) is the perturbation amplitude, 0p is the initial pre- 
cession phase, Tp is the precession period, cOp - 2n/Tp and 
kp = LOp/VA- The time t is measured in the rest frame of the 
beam. The tilde in x and ^ indicates that the orbital motion 
has not yet been taken into account (here the axes x and ^ are 
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assumed to be parallel to the axes x and ^, respectively; see 
Fig. Obviously, ^ - ^, and the tilde can be omitted in all 
expression with ^. The R(^) can be chosen to be 



(5) 



Here, a - /?o/(2tanQp), and R() is the maximum amplitude of 
the perturbation. 

The form of the function ^(f) depends on the evolution of 
the speed j6c of the plasma condensation injected into the per- 
turbed beam. From Eq. 0, the components of fi^ are: 



ySf = P^idRIA^ cos(wpf - fep^ + 0p)- 

- R{^) (w - kp^) sin(wpf - k^^ + 0p) 

Pi = P^midi^) sin(wpr - k^^ + 0p)+ (6) 
+ R{^) (w - kp^) cos(wpf - k^^ + 0p) 
= d^/df. 



so that 



,6? = (d^/df)'Ci + R\^) (Wp - k^d^lAtf , (7) 



whereCi = (d7;/d^)^-i- 1, and remembering that -PI+PI + 
/??. Equation is a quadratic equation with respect to d^/df, 
and it can be written in the form of 



A2(d^/df)^ + Ai(d^/df) + Ao = , 



(8) 



with A2 = 1 + R'-{^)kl + (d«/d^)2, Aj = -lR^(^)io^k^ and 
Ao = 7?^(^(;)Wp - y6^. The component speed jSc can be con- 
veniently expressed as a function of ^ or ^(f) (it can also be 
viewed constant, in the simplest case). The jet side solution of 
Eq. Q can be obtained for the time f(^) in the rest frame of the 
beam, under the condition of initial acceleration of the beam 
(d^/c/f > 0, for^ — > 0+). The formal solution of Eq. (|8} is 



df 



1 
-I- 



■{R^(^)u,^k^^ 



dR. 



R\Ok'^pl + (PI- R\0<^1){1 + ^) 



1/2 



(9) 



and t(^) is then given by 



m = fo + 



kpOjpRH^) + [AlPl - ulR\^)Ci]'l^ 



d^. (10) 



In general, both Wp and k^ can also vary. Their variations should 
then be given with respect to ^, in order to preserve the solution 
given by Eq. J10> . Numerical methods should be used other- 
wise. The solution determined by Eq. JlOl i can be found for all 
R{^) not exceeding the limit of 



^(^)<y8c[w^-fe^?cosV(^)] 



-1/2 



(11) 



The beam pitch angle i/'(^) = arctan(d/?/d^) in Eq. M \\ is 
measured "locally" (as opposed to the "global" pitch angle 
i^g(^) = arctan[7?(^c)/f]), and it describes the evolution of the 
amplitude of the beam perturbation. For a jet with constant 



beam pitch angle (i/'(^) - ^ - const), the inequality [H] can 
be rewritten as: 

-1/2 



^<Pc [(j^l tan^ (A - kjpl sin^ ij/]' 



(12) 



The speed p^ and the local pitch angle t// should satisfy the con 
dition: 

Va 



Pc cos tp < 



(Vl + 1)« 



(13) 



If inequalities II 1H13I are satisfied, the solution provided by 
Eq. ilOi can be used to represent the form of the function ^(f) 
and describing the kinematics of the beam. 

2.4. The binary system 

The trajectory of a superluminal feature travelling in the per- 
turbed beam is modified also by the orbital motion of the two 
black holes in the binary system. This motion can be taken 
into account by solving the restricted three-body problem for 
the BBH-AD system (LaskarT^ Laskar & Robutel[T995}. 
The orbital plane is {xOO, and the coordinate system is cen- 
tered at the mass center of the system. The orbit is given by 
r - p/[l + e cos^(f)], where (f is the true anomaly, p is the 
semi-latus rectum and e is the numerical eccentricity of the or- 
bit. The resulting coordinates (x\, ^\) of the black hole Mi are 
given by 



Xiit) = 

m 



M2 



M\ + M2 I + e cos((p(t)) 
M2 p 



cos(<p(t)) 



sm(<p(t)) 



(14) 



Ml + M2 I + e cos((^(f)) 

The beam trajectory can be reconstructed from Eqs. (0} and 
(I14> . which gives 



X=X+Xl(t) COS(dL»of - ^o^(f) + 4>o) 

^ = i + ^i(t) sinicoj - ko^(t) + 



(15) 



where 0o is the initial orbital phase, = ^n/To, ko = oJo/Va, 
and To is the orbital period. This system of equations describes 
completely the kinematics of the perturbed beam. Similarly to 
the solution obtained for Eq. ^(f) can now be determined 
by forming pl from the coordinate offsets given by Eq. (I15> and 
solving for d^/df, so that (for the jet side solution) 



(Ai2-4AiAo)i/2-A 



■d^ 



(16) 



The coefficients Aq, Ai, and A2 are given explicitly in 
Appendix II in Britzen et al. (^001 1. 

2.5. Properties of the t(^) solution 

The strongest effect on the form of the beam trajectory is 
produced by three model parameters: the beam speed p^, the 
Alfven speed Va, and the precession period Tp. The orbital pe- 
riod in the BBH system affects only the position of the foot- 
point of the jet. Since VLBI measurements are made with re- 
spect to the footpoint position. To would not have an effect on 
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fitting the observed trajectories of the jet components (the or- 
bital period is constrained by the variability of the optical emis- 
sion). It is easy to see that the remaining parameters (/?o, fip, 
I and 6) describe only the orientation and shape of the overall 
envelope within which the perturbed beam is evolving. The t(^) 
solutions described by Eqs. ilQ\ and^]can degenerate or be- 
come non-unique with respect to y6c, Va, Tp, if these parameters 
are correlated. The combinations of j8c, Tp and Va describing 
the same beam trajectory must not depend on the rest frame 
time t or beam location ^. However, the solution given by equa- 
tion[lO|implies that the beam trajectory would remain the same 

Obviously, Eq. (^} depends on both t and ^. The dependence 
vanishes only when t — > ^/c. The precession period then be- 
comes a complex number, and the resulting combination of the 
parameters does not describe a physically meaningful situation. 
This implies that the model does not become degenerate or 
non-unique with respect to j6c, Tp and Va, and each physically 
possible combination of these parameters describes a unique 
beam trajectory. 

2.6. Specific properties of the model 

The model introduces two new properties in the description of 
the BBH-AD system: 1) the Lorentz factor of the relativistic 
beam is time variable; 2) the initial perturbation of the mag- 
netic lines in the beam dissipates on a timescale fbeam- These 
modifications are introduced in order to accommodate for the 
kinematic properties of the jet in 3C 345 in which the curvature 
of component trajectories decreases at larger distances, and ac- 
celerated motions are strongly implied by the observed trajec- 
tories and apparent speeds of several jet components (LZ99l. 
The beam accelerates from Fmin to r^ax over a rest frame time 
interval face, and 



rc(f) = 



(t/kcc) 



1 + f/fa. 



(18) 



The perturbation dissipates at a characteristic time fbeam- 
The amplitude of the perturbation given by equation Q is mod- 
ified by an attenuation term (characteristic damping time of the 
perturbation), so that 



(19) 



The model takes account of an arbitrary source orientation, 
by introducing a position angle, 0^ of the ^-axis in the picture 
plane of the observer. 

2. 7. Observed quantities 

Corrections for the relativistic motion of the beam must be 
made, in order to enable calculations of physical quantities in 
the observer's frame. This requires a knowledge of the angle 
6{t) between the instantaneous velocity vector of the compo- 
nent and the line of sight (LOS). Since the observed angular 



size of the BBH orbit is small (<k 1 milliarcsecond), the LOS 
can be assumed to lie in the plane parallel to (xO^), making an 
angle i with the ^ axis. In this case. 



/dy df \ 1 

cos 6(t) - \ — sm i H cos l — 

\df df 1/3, 



(20) 



in the observer's frame (Camenzind & Krockenberger [T992t . 
with the resulting Doppler beaming factor 6^(1) - {r^.[l - 
ySc cos 0( ?)])"' The observed flux density of the moving com- 
ponent is 



1 



5e = -2 4(0'"""(l+z) 

Df 



jcdV, 



(21) 



where z and Dl are the redshift and luminosity distance of the 
source, is the emissivity of the component and el's is the 
synchrotron spectral index. The shortening of time in the ob- 
server's frame is given by 



fobs = (1 +Z) 



X re( 



df 



c(f)<Jc(f) ' 

and the component travels over the distance 



^trav(f) = (1 + Z)^ 



3. The quasar 3C345 



ySc COS ^(f) 



-df . 



(22) 



(23) 



The 16™ quasar 3C345 exhibits remarkable structural and 
emission variability on parsec scales around a compact un- 
resolved radio core. The total radio flux density of 3C345 
has been monitored at 5, 8, and 15 GHz (Aller, Aller & 
Hughes I2003J . and at 22 and 37 GHz (Teriisranta et al. lT^ . 
The source has also been monitored with the Green Bank 



Interferometer at 2.7 and 8.1 GHz (Waltman et al. [19911 . The 
observed variability of the optical emission is possibly quasi- 
periodic with a period of »;1560 days (Babadzhanyants & 
Belokon' 1984; Kidger 1989), although it has been suggested 
that the light curve may originate from a non-linear and non- 
stationary stochastic process (Vio et al. I1991> . A strong opti- 
cal flare observed in 3C 345 in 1991-92 (Schramm et ai. \TW^ 
Babadzhanyants, Belokon' & Gamm 1995J was connected 
with the appearance of a new superluminal jet component, C7 
(|L96 ). LZ99 report possible 3.5^ year quasi-periodicity man- 
ifested by spectral flares of the radio emission and appear- 
ances of new jet components. Historical optical and 8 GHz ra- 
dio lightcurves are shown in Fig.|3]together with the epochs of 
the spectral flares and component ejections. 

The evolution of the parsec-scale jet of 3C 345 has been 
studied extensively with the VLBI: see Unwin et al. (I1983> ; 
Biretta, Moore & Cohen (Il986i) ; Zensus, Cohen & Unwin 
(1995, hereafter IZCU95t ; 06l Ros, Zensus & Lobanov 
(2000, hereafter lOO^; Klare ("2003). The relativistic jet model 
(Blandford & Rees 1978j has been applied to explain the emis- 
sion and kinematic changes observed in the parsec-scale struc- 
tures. Using the X-ray data to constrain the jet kinematics, 
IZCU95I and Unwin et al. (1994) derive physical conditions 
in the jet from a model that combines the inhomogeneous jet 
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Fig. 3. Optical and radio variability of 3C 345. Shown are the epochs of spectral flares (ILZ99> that have resulted in appearances 
of new superluminal components in the jet. The radio data are from the Michigan monitoring program ( Aller, AUer & Hughes 



EOnS*. The optical data are from Kinman et al. (il968< : Smyth & Wolstencroft Q9T0i; Lii C1972.) : McGimsey et al. (IT975t 
Pollock et al. iWMy Angione et al. ([1981); Kidger ( fl^SSt Webb et al. (IT988t : Kidger & de Di ego (l990l : Vio et al ST99i\ : 
Schramm et al. f 1993 1; Babadzhanyants, Belokon' & Gamm J1995> : Belokon' & Babadzhanyants J1999> . All optical data have 
been converted into the B magnitude scale. 



model of Konigl ( I1981> for the core with homogeneous syn- 
chrotron spheres for the moving jet components (Cohen 1985 1. 
The evolution of the core flux density is well represented by a 
sequence of flare-type events developing in a partially opaque, 
quasi-steady jet (LZ99l. Spectral properties of the parsec-scale 
emission have been studied in several works CL96i, Lobanov, 
Carrara, & Zensus lT997l Lobanov [T998bl . Unwin et al. {1997 1 
have analyzed a correlation between the X-ray variability and 
parsec-scale radio structure of 3C 345. 

3. 1 . Compact jet in 3C 345 

Figure 0] shows VLBI images of the jet in 3C 345 at 5 and 
15 GHz, with the trajectories of eight superluminal components 
identified and monitored in the jet in 1979-1999. The morphol- 
ogy of 3C 345 is typical for a core dominated source: a bright 
and compact core responsible for 50-70% of the total emis- 
sion, and a compact jet that contains several enhanced emis- 



sion regions (jet components) moving at apparent speeds of up 

to ySapp ~ 20/1-1. 

The compact jet in 3C 345 extends over ~ 25 mas, which 
corresponds to a projected distance of ^95 (/z sin^jet)"' pc 
(a;800/!-' pc, assuming the most likely viewing angle of the 
jet fljet = 7°) . The jet has two twists at the separations of ail 
and x5 mas. The twists appear to retain their positions in the jet 
over the entire duration of the VLBI monitoring. Such stability 
of the twists suggests that they are likely to reflect the three- 
dimensional structure of the jet. The unchanging positions of 
the twists are also consistent with the similarities of the ap- 
parent trajectories of the jet components at separations >1 mas 
(Fig-0- The two-dimensional trajectories observed in several 
jet components follow similar, overlapping tracks, and outline 
a continuous picture of a curved jet. The separation of 20 mas 
from the core can mark a transition point from compact jet to 
large-scale jet. At this separation, the jet finally turns to the 
direction (p » -40° observed in the large-scale jet of 3C 345 
(Kolgaard, Wardle & Roberts 1989) . 
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Fig. 4. Parsec-scale jet in 3C 345 at 5 and 15 GHz. Contours are drawn at -1,1, '\/2,2,... of the lowest contour (2.5 mJy/beam, in 
both images). Curved lines mark the trajectories of superluminal features C 1-C7 identified in the jet. The component trajectories 
are polynomial fits to the position measurements from the 1979-1998 VLBI monitoring data (iiZCU95. 1^96. .ROQj . The positions 
of C7 measured at 22 GHz are shown in the inset. 



The measured positions of the component C7 (Fig.0} pro- 
vide an accurate account of the jet kinematics within ~2mas 
distance from the core. The large amount and good quality 
of the data at 22 GHz are sufficient for determining the com- 
ponent trajectory without using measurements made at lower 
frequencies. The total of 18 VLBI observations of 3C345 at 
22 GHz are combined for the epochs between 1989 and 1998 
(L96 ROO Klare 2003 1. For several of these observations, new 
model fits have been made for the purpose of obtaining consis- 
tent error estimates for the positions and flux densities of C7. 
Table|3]presents the resulting database of the positions and flux 
densities of C7 at 22 GHz. 

C7 was first rehably detected in 1991.5 at 22 GHz at a 
core separation of ©1991.5 ~ 100 //as {L96j Reports of ear- 
Uer detections in 1989.2 at 100 GHz (Baath et al.[T992j and 



22 GHz (!L96^ are inconclusive because of the likely confusion 
and blending with the nearby component C6 and the core itself. 



Following the procedure established in IZCU951 and lL96l 
polynomial functions are fitted to the component coordinate 
offsets x{t) and y(t) from the core position which is fixed at the 
point (0,0). From the polynomial fits, the proper motion fi{t) = 
[(dx/df)^+(dy/df)^]'^^, apparent speedySapp(0 = /"(0£*l(1+z)"' 
and apparent distance &app(t) - J i^(t')dt' travelled by the 
component along its two-dimensional path can be recovered. 
For the purpose of comparison to the previous work of L96 and 
Lobanov & Zensus ( 1996 1, the detection reported for the epoch 
1989.25 is included in the polynomial fitting, but all results for 
the epochs before 1991 should be viewed as only tentative. 
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Fig. 5. Kinematic properties of the jet determined from the observed trajectory of C7. In all panels, the solid line represents a 
model with Fc? increasing linearly with time from Fc? = 2.2 in 1990 to Fc? = 10.6 in 1998. The lower boundary of the shaded 
area corresponds to Fc? = F^in = (filpp + 1)'^^, the upper boundary is drawn for Fc? = 11. Arrows indicate the direction of 
increasing Fcv from Fmin to 1 1 . Individual panels show: (a) Measured apparent speed ySapp (dashed line) and evolution of the 
Lorentz factor Fj. (b) Angle between the velocity vector of C7 and the line of sight, (c) Doppler factor, (d) Distance travelled by 
C7 in the rest frame of the jet. The derived T^trav > jScvAfobs because of the time contraction in the observer's frame of reference. 



3.2. Kinematics of tlie jet component C7 

Kinematic properties of C7 are presented in Fig. |5l The ap- 
parent speed jSapp of C7 increased from x 2 c in 1989 to 
j6max ~ 13 c in 1997 (the decrease of ySapp after 1997 is an ar- 
tifact of the polynomial fit). The component reached the min- 
imum apparent speed of ySmin ~ 0.5 c in 1990.5. Such large 
changes of jSapp indicate that C7 has moved along a substan- 
tially curved three-dimensional path, with the bulk Lorentz fac- 
tor most likely changing with time. If the Lorentz factor is con- 
stant, then Fc7 > 13 (to satisfy the observed jSapp)- If rmathrmC? 
varies, the lower boundary of the variations is given by the min- 
imum Lorentz factor F^in = Wlpp(t) + 1]'^^. which describes 
the motion with the least kinetic power. In Fig.|5j F^ax = 15.5 
is set, and the kinematic parameters of C7 are described for 
all allowed variations of Fc7(f) (Finin(f) < Fc7(f) < F,nax)- As 
seen in Fig. [S] neither of the two boundaries provides a sat- 
isfactory description of the component kinematics. The case 
Tc? = Tmax results in uncomfortably small Ociit) and dispro- 



portionally large travelled distance /?trav(0- With Fc? - ^mm(t), 
the resulting variations of Oci exceed 50°, which is difficult 
to explain. A plausible evolution of the rest frame speed of 
C7 can be described by a linear increase of the Lorentz fac- 
tor from Fc7 = 2.2 in 1989 to Fcv = 15 in 1998. In this 
case, none of the kinematic parameters change drastically in 
the course of the component evolution: Oqi varies within x 10°, 
and /?tiav ^ lOpc. This indicates that the motion of C7 is likely 
to be accelerated during the initial stages of its evolution, at 
distances of < 2 mas from the core. At such distances, acceler- 
ation of the component rest frame speeds has also been reported 
in several other jet components in 3C 345 (Lobanov & Zensus 
II996IILZ99> . Consequently, this must be a general property of 
the jet plasma on these spatial scales, and it should be taken 
into account in the modelling of the jet kinematics. 
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Table 1. Positions and flux densities of C7 at 22 GHz L 



fobs 


S 




Ax 


Ay 




[year] 


[Jy] 


[mas] 


[mas] 


1989.25^ 


0.19 ± 


0.02 


0.04 


±0.02 


0.01 ± 


0.01 


1991.46 


0.80 ± 


0.04 


0.09 


±0.05 


-0.01 ± 


0.01 


1991.86 


3.29 ± 


0.49 


0.12 


±0.01 


0.03 ± 


0.01 


1992.45 


6.20 ± 


0.36 


0.18 


±0.01 


0.00 ± 


0.01 


1992.86 


5.49 ± 


0.30 


0.28 


±0.01 


0.01 ± 


0.01 


1993.13 


3.60 ± 


0.10 


0.34 


±0.01 


-0.02 ± 


0.01 


1993.72 


2.23 ± 


0.12 


0.60 


±0.04 


-0.10± 


0.02 


1 QQA A ^ 


1.04 ± 


0.10 


0.90 


±0.09 


-0.12 ± 


0.08 


1995.85 


1.21 ± 


0.06 


1.20 


±0.08 


-0.10± 


0.08 


1996.42 


0.88 ± 


0.05 


1.36 


±0.08 


-0.07 ± 


0.07 


1996.82 


0.69 ± 


0.04 


1.55 


±0.08 


-0.08 ± 


0.12 


1997.40 


0.64 ± 


0.10 


1.87 


±0.10 


-0.22 ± 


0.13 


1997.62 


0.59 ± 


0.07 


2.02 


±0.10 


-0.24 ± 


0.12 


1997.90 


0.47 ± 


0.07 


2.10 


±0.10 


-0.30 ± 


0.11 


1998.02 


0.47 ± 


0.08 


2.02 


±0.25 


-0.24 ± 


0.14 


1998.99 


0.50 ± 


0.10 


2.63 


±0.20 


-0.15 ± 


0.14 


1999.12 


0.28 ± 


0.10 


2.69 


±0.20 


-0.36 ± 


0.14 


1999.64 


0.17 ± 


0.20 


2.87 


±0.20 


-0.34 ± 


0.15 



Note: t - detection is likely to be confused due to blending between 
the core and the preceeding jet component C6. References: 1989.25- 
1993.72 i L96j : 1994.45 (Leppanen, Zensus & Diamond .1995 ): 
1996.42-1997.40 jM^ : 1997.62-1999.64 (Klare l^TO . 



3.3. Long term trends in the jet: evidence for 
precession 

The observed consistency of the two-dimensional component 
tracks plotted in Fig. 0] suggests that all jet components in 
3C345 travel along similar trajectories. There is however a 
substantial amount of evidence for global evolution of the en- 
tire jet. In Fig.|6j the apparent accelerations, dyu/d0app, of the 
components are plotted against the respective component ori- 
gin epochs taken from L96 The apparent accelerations show 
a steady increase, following the order of component succes- 
sion, which is a strong indication for long-term changes in the 
jet. While it is conceivable that the younger components are 
intrinsically faster, this would require the presence of a non- 
stationary process affecting significantly the particle accelera- 
tion in the nuclear regions of 3C345. The more likely expla- 
nation of the observed long-term trend of the component ac- 
celerations is precession of the jet axis. In this case, variations 
of the component ejection angle alone are sufficient to explain 
the observed trend: the components ejected at angles closer 
to flfast - 1/rb would have higher dyu/d0app. In the course of 
the precession, the ejection angle would change steadily, and 
the measured dyu/d©app should show sinusoidal variations. The 
rather large errors in the measured dyu/d©app in Fig. |6] allow 
only a very rough estimate of the precession period Tpiec, which 
should be in the range of 850-3670 years, to be made. 

Periodic changes are also evident in the variations of the jet 
position angle measured at 22 GHz for successive jet compo- 
nents as they reach a ^0.5 mas separation from the core. The 
resulting position angles are plotted in Fig.0 The position an- 
gle varies with a period of Pshort ~9.5 years and an amplitude 
of Ashoit ~20°. A longer term trend is suggested by the slope 
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Fig. 6. Apparent acceleration d///d0app of the component 
proper motions, plotted versus the component epochs of ori- 
gin. The origin epochs are obtained by back-extrapolating the 
observed trajectories of the jet components (L96). The shaded 
area cover the range of precession periods which would repro- 
duce the observed trend of dyL(/d0app. 



in the linear regression of the position angle data, which im- 
plies a rate of change in the position angle of =s0.4°/year. The 
time coverage of the data is too short to provide an accurate 
estimate of the period of the long term variations. A formal 
fit constrains the period to within a range of 300-3000 years, 
with a mean value of Piong ~ 1600 years and an amplitude of 
Along ~ 10°. The shaded area in Fig.0represents a combined fit 
to the data by the two periods. The estimate of Piong obtained 
from the change in position angle falls within the range of pe- 
riods implied by dju/d0app, and it can therefore be associated 
with the precession. The nature of the short-term changes of the 
jet position is unclear, but they could be related to the rotation 
of the accretion disc or the jet itself or to the orbital motion in 
the BBH system. 

4. Application of the BBH model 

The BBH model applies a single model framework to explain 
the emission and kinematic properties of the jet in the opti- 
cal and radio regimes. The evolution of superluminal features 
observed in the radio jet is explained by a precession of the 
accretion disk around the primary black hole Mi, and the char- 
acteristic shape of the optical light curve is produced by the 
orbital motion of the black hole Mi . 

The different characteristic timescales and amplitudes of 
changes in the radio and optical regimes imply a complex 
composition of the emitting material. The optical emission of 
3C 345 during the strong flare in 1990-93 varied by » 2.5 mag, 
with a typical timescale of ~0.6yr for individual sub-flare 
events (see Fig. I15> . The radio emission of the jet compo- 
nent C7 associated with that flare showed only a single rise 
and decay event, changing by a factor of ~ 15 on a timescale 
of ~6 years. This difference in the behavior of the optical and 
radio variability can be reconciled by assuming that the rel- 
ativistic plasma responsible for the optical emission is in- 
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Fig. 7. Position angle of different jet components measured at 
22 GHz at 0.5 mas separation from the core. Short-term vari- 
ations, with a period of ~9.5 years, are evident. A linear re- 
gression (dotted line) suggests the presence of a long-term 
trend, with a annual change rate of 0.4+0. 3°/year. The shaded 
area shows the best fit by the two periodical processes com- 
bined, with the estimated periods Pshon - 9.5 + 0.1 years and 
^long = 1600 + ISOOyears. The solid line represents the trend 
in the position angle due to the long-term periodicity. 



jected into the beam during a time Aropt which is much shorter 
than the injection timescale Arjad for the plasma responsible 
for the radio emission. For C7, the plasma responsible for the 
optical emission is modelled by a point-like component and 
the plasma responsible for the radio emission is modelled by 
a temporally extended component, so that Axopt < O.OlArrad- 
In the BBH frame, the radio component emits over a period 
of =s 1.9 yr, with component ejections occurring approximately 
every 4 years. A general sketch of the composition of a flare is 
shown in Fig.|8l 

The application of the model consists of two parts. In the 
first part, the observed trajectory and flux density of the jet 
component C7 are reconstructed, using the precession of the 
accretion disk around the black hole Mi, and solving formally 
the problem defined by Eq. (|8}. In the second part, the optical 
light curve of 3C 345 is calculated from the orbital properties 
of the BBH system provided by the solution of Eq. il6\ . 



4.1. Precession and geometry 

Basic geometrical parameters of the system are recovered by 
fitting the precession model to the observed two-dimensional 
path of C7. These include 
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Fig. 8. Composition of a nuclear flare. The flare begins with 
an injection of highly energetic particles emitting optical syn- 
chrotron radiation. The timescale of the injection Aropt is com- 
parable to the synchrotron loss time of the optically emitting 
plasma. Rapid energy losses and injection of radio emitting 
plasma on the timescale Ar^d lead to progressively decreas- 
ing mean energy of the plasma and subsequent propagation of 
the flare to radio wavelengths. The optically emitting plasma 
is represented in the model by a point-like component injected 
into the beam, while the radio emitting plasma is temporally 
resolved, with the resulting Aropt < O.OlAriad 

Va of the magnetic field perturbation. This speed is determined 
by fitting simultaneously the geometry of the component path 
and the orbital motion in the binary system. 

4.1 .1 . Velocity of the jet component 

The velocity of the jet component is constrained by fitting the 
coordinate off'sets, x(t) and y(f), of the component from the po- 
sition of the radio core. In 3C 345, the apparent motion of all 
jet components is non-linear and shows accelerations and de- 
celerations JZCU95I IL961 ILZ99t . To account for this, a vari- 
able Lorentz factor must be introduced, using Eq. J18I I. With a 
variable Lorentz factor, the x(t) and y{t) off'sets are fitted simul- 
taneously, and a three-dimensional trajectory of the component 
is recovered. This procedure yields the following parameters: 

rmin - initial Lorentz factor of the component; 
Tmax - final Lorentz factor of the component, after the accel- 
eration; 

face - characteristic acceleration time. 



I - viewing angle of the jet axis (Fig.O; 

"Pproj - position angle of the jet axis in the plane of the sky; 4.2. F/ux cfens/fy changes in the jet component 



Tp - precession period; 
Qp - opening half-angle of the precession cone; 
(^p - initial phase of the precession; 
Ro - maximum amplitude of the perturbation; 
fbeam - characteristic dampening time of the perturbed beam. 

It should be stressed that none of the parameters listed 
above depend on the velocity of the moving component. 
However, the shape of the path depends on the Alfven speed 



In addition to the geometrical and kinematic factors, the flux 
density evolution in the optical and radio regimes is charac- 
terized by the synchrotron energy loss timescales Aropt and 
Ariad and by the duration Arthick during which the relativistic 
plasma responsible for the flare remains optically thick in the 
radio regime. The radio lightcurve of C7 constrains uniquely 
the viewing angle of the jet. The information contained in the 
radio lightcurve is fundamental for determining the duration of 
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the radio emitting stage, At^-^^, of the flare and ensuring the 
consistency between the precession fit and the BBH fit. 

4.3. BBH system 

The main effect of the BBH system is the orbital motion of 
both black holes around the center of mass of the system. This 
motion explains the observed optical light curve, the fine struc- 
ture of the parsec-scale jet, and peculiarities of component mo- 
tions around the general trajectory determined by the preces- 
sion. The fit by the BBH system yields the following parame- 
ters: 

Ml - mass of the primary black hole containing the accretion 
disk; 

M2 - mass of the secondary black hole; 

To - orbital period in the BBH system; 

e - eccentricity of the orbit; 

to - injection epoch of a new jet component; 

4>o - initial orbital phase at the injection epoch. 

4.4. Consistency between the precession and BBH fits 

The BBH model, which is obtained by fitting a point-like com- 
ponent to reproduce the observed optical lightcurve, must be 
consistent with the results derived from the precession fit. In 
other words, the evolution of the extended, radio emitting, com- 
ponent determined from fitting the BBH model must reproduce 
the trajectory and kinematic properties of the radio emitting 
plasma obtained from the precession fit. 

Consistency between the fits by the precession and orbital 
motion is ensured under the following five conditions which 
connect the formal fitting with the specific properties of the 
two-fluid model: 

1 . The characteristic Lorentz factor obtained from fitting the 
optical light curve by the BBH model must be larger than 
the corresponding Lorentz factor obtained from fitting the 
trajectory of C7 by the precession model. This is simply due 
to the fact that the rest frame trajectory in the BBH model 
is longer than the trajectory in the precession model (see 
Fig. ^5 implying that the component has to move faster 
in the BBH model to reproduce the same observed two- 
dimensional trajectory. 

2. The optical component is assumed to be temporarily unre- 
solved, but not the radio component. This condition is re- 
quired to explain the differences in the emission variations 
observed in the optical and radio regimes. 

3. The ejection epoch of the optical component precedes that 
of the radio component. This condition reflects the nuclear 
flare structure assumed (Fig.|8}. 

4. For a given mass Mi of the primary black hole, the combi- 
nation of the mass of the secondary M2, orbital period To 
and acceleration time face must reconcile satisfactorily the 
rest frame time derived from the precession model with the 
rest frame time obtained from the BBH model (see Fig.O. 

5. A single value of the Alfven speed must reproduce the ob- 
served shape of the component trajectory obtained from the 



precession fit and satisfy the orbital period obtained from 
the BBH fit. 

These five conditions ensure that the kinematic and flux 
density variations of the radio component recovered from the 
BBH and precession fits are consistent, and thus the BBH sce- 
nario explains the radio and optical components simultane- 
ously. 

5. Properties of the BBH system in 3C345 

The best fit parameters of the precessing jet and the BBH sys- 
tem in 3C 345 determined using the model developed in Sect.|2l 
and the fitting method described in Sect. ^ are presented in 
Table |2] The procedural steps through which the fit was ob- 
tained are summarized below. 

5.1. Geometrical and precession properties 

The position angle of the jet axis in the plane of the sky is 
determined directly from the VLBl images of 3C 345, which 
yields 

^proj — 95.5 . 

The injection epoch to does not have to coincide with the 
epoch of first detection of C7, since the temporal and linear 
scales for plasma acceleration may be comparable to, or even 
larger than, the smallest angular scale at which the component 
becomes detectable in VLBI images. The most accurate esti- 
mate of the epoch of origin of C7 is provided by the onset of 
the radio flux increase (see Fig.|3}, which yields 

to = 1991.05 

This value indeed provides the best fit for both the optical 
lightcurve and the evolution of C7. 

The jet viewing angle, l, opening half-angle of the preces- 
sion cone Qp, and initial precession phase, <pp, are obtained by 
fitting the observed two-dimensional path of C7. This fit con- 
strains uniquely the initial precession phase at to 

0p = 100° . 

For fitting the jet viewing angle, an initial guess is pro- 
vided by the best fit values (t - 5-8° obtained from mod- 
elling the kinematic JZCU95t . spectral (|LZ99l and opacity 
(Lobanov 1 1998 a^ properties of the compact jet in 3C345. A 
family of solutions across the t = 5-8° interval is presented 
in Table|3] together with corresponding T^ and fbeam calculated 
for Va = 0.078 c. (note that determination of t and Qp does not 
depend on the knowledge of Va)- The best fit to the observed 
path of C7 is provided by 

i = 6.5°, Qp = 1.45°. 

At this stage, the maximum perturbation amplitude, Rq, and 
the timescale, fbeam of damping the perturbed beam can be con- 
strained. An initial guess for/?o (0.5-0.8 pc) is provided by the 
largest deviation of the positions of C7 from a straight line de- 
fined by Tproj . The best fit for the observed path of C7 is ob- 
tained with 

Ro^l.Opc, fbeam = 180yr. 
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Table 2. Precession, orbital and flare properties of the BBH system in 3C 345 



Geometrical parameters 



Viewing angle of the jet, 


L 


6.5° 


Position angle of the jet in the plane of the sky, 


i pi-oj 


-95° 


Opening half-angle of precession cone. 




1.45° 


Precession period of accretion disk. 


T 


2570 yr 


Characteristic dampening (comoving) timescale of perturbation, 




180yr 


Maximum perturbation amplitude. 




1 pc 


Kinematic parameters 


Alfven speed. 


Va 


0.078 c 


Initial Lorentz factor of radio emitting plasma. 


r 

^ mm,r 


3.5 


Initial Lorentz factor of optically emitting plasma. 


r 

^ mm,o 


12 


Terminal Lorentz factor of jet plasma. 


r 

A max 


20 


Plasma acceleration (comoving) timescale. 


face 


165 yr 


Parameters of the binary system 


Primary black hole mass. 


Ml 


7.1 X 10** Mq 


Secondary black hole mass, 


M2 


7.1 X lO'* Mq 


Orbital period. 


T 


480 yr 


Orbital separation. 




0.33 pc 


Orbital eccentricity. 


e 


<0.1 


Characteristic rotational period of accretion disk. 




240 yr 


Flare parameters and initial conditions 


Plasma injection epoch, 


Anj 


1991.05 


Initial precession phase. 


0p 


100° 


Initial orbital phase. 


<Po 


125° 



Table 3. Precession solutions for ( = 5-8° 



Parameter 


Value 




5 


6 


6.5 


7 


8 




1.3 


1.4 


1.45 


1.5 


1.6 


'beam [y] 


220 


200 


180 


165 


150 


Tp [yr] 


3320 


2780 


2570 


2370 


2070 



Note: Bold face denotes the best fit solution. Tp is calculated with the 
best fit value of the Alfven speed Va = 0.078 c. 



5.2. Kinematic fit 

After the two-dimensional path of C7 has been fitted by the ge- 
ometrical and precession parameters, the precession model is 
applied to fit the trajectory and flux density evolution of C7, 
which yields r^in - 3.5, r^ax = 20 and face - 165 yr These pa- 
rameters are independent from the choice of Va- It should also 
be noted that the Lorentz factors derived above are related to 
the radio emitting plasma. The optically emitting plasma may 
have a different range of Lorentz factors. 



Table 4. Precession solutions for Va - 0.05 c-0.15 c 



Parameter 


Value 


Va [C] 


0.050 


0.075 


0.078 0.100 


0.125 


0.150 


T, [yr] 


4060 


2690 


2570 1970 


1510 


1230 



The fit for Tp depends on the Alfven speed Va in the jet. For 
each value Va an appropriate Tp can be found from the preces- 
sion fit that would reproduce exactly the same two-dimensional 
trajectory. This family of solutions is given in Table ^ The 
value of Va is constrained after a similar family of solutions 
is obtained for the BBH fit to the optical light curve. A unique 
solution for Va is then found by reconciling the values of Tp 
and To obtained from the BBH and precession fits. 



5.3. Consistency between tlie radio and optical fits 

In order to ensure consistency between the precession and BBH 
fits, and to bootstrap together the fit of the trajectory and kine- 
matics of C7 and the fit of the optical light curve, two gen- 
eral model classes are considered for describing the proper- 
ties of the radio and optically emitting plasmas, in terms of the 
minimum (or initial) Lorentz factors of the two plasmas (these 
scenarios are based on the consistency conditions described in 
Sect.iS: 

Model 1: Fmin.rad < A min, opt 5 
Model 2: Fmin.rad - rmin.opt- 

In both cases, the radio emitting plasma is injected immediately 
after the optically emitting plasma (as illustrated in Fig.|8|l. 



For Model 1,F, 



mill, rad 



4 and F, 



mill, opt 



12 are chosen from 



the precession fit (note that this corresponds to injection speeds 
of Vmin^rad - 0.9682 c and Vmin.opt = 0.9965 c). For Model 2, 
rmin.rad = rmin.opt = 12 is assumcd. Both sccnarios fit well 
the optical light curve, but Model 1 reproduces better the kine- 
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Fig. 9. Rest frame time recovered from the precession fit (de- 
scribing the trajectory of the VLBI component) and the BBH 
Models 1 and 2 (describing the optical variability) under the 
consistency constraints discussed in Sect. 14.41 The injection 
epoch finj - 1991.05 (for which fiest = is required by the 
precession fit) is reproduced by the BBH Model 1 and cannot 
be recovered from the BBH Model 2. The agreement between 
the precession fit and Model 1 ensures that the BBH scenario 
provides a single framework for describing simultaneously the 
observed trajectory and flux density changes of C7 and the ob- 
served variability of the optical emission of 3C 345 

matics of C7, particularly at the beginning of the component 
evolution. 

Fig-Elcompares the rest frame time obtained from the pre- 
cession fit and the BBH fit by the two models. The injection 
epoch finj = 1991.05 obtained from the precession fit is recov- 
ered by Model 1, but cannot be reproduced by Model 2. 

For the earliest epochs of the evolution of C7, Model 2 does 
not reproduce satisfactorily the separations of C7 from the core 
(Fig. I10> and the flux density changes (Fig. llH . The discrep- 
ancy between the fits by the precession model and Model 2 
is particularly visible in the speed evolution obtained from the 
precession and BBH fits (Fig.ll2>. Model 1 is therefore a better 
counterpart to the precession fit described in Sect. 15. II 

Figure [O] compares the precession and the BBH Model 1 
fits to the observed two-dimensional path of C7. The agree- 
ment between the fits is excellent. The two fits difi'er by less 
than 0.02 mas (Fig.ll4>. This deviation is caused by the orbital 
motion of the black hole ejecting the jet, which is not accounted 
for in the precession fit. 

5.4. Optical light curve and orbital period 

Fitting the optical light curve constrains the orbital period To, 
orbital eccentricity and initial orbital phase <po in the BBH sys- 
tem, under the conditions of consistency described in Sect. 15. 31 
The fit to the optical curve requires a nearly circular orbit 
(e < 0.1), constrains uniquely the initial orbital phase 

<Po = 125° 

and produces a family of solutions for To depending on the 
value of Va- These solutions are presented in Table |5] for the 
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Precession modei 
BBH modei 1 
BBH modei 2 
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1994 1996 

Observed Time 



Fig. 10. Evolution of the separation of the jet component C7 
from the core. The fit by the precession model is consistent 
with the fit by BBH Model 1. BBH Model 2 fails to reproduce 
the earliest position measurements of C7. 
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Fig. 11. Radio flux density of C7 at 22 GHz. The fits by the pre- 
cession model and BBH Model 1 are consistent. BBH Model 2 
cannot reproduce the onset of the radio emission. 
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Fig. 12. Apparent speed evolution of the jet component C7 re- 
covered from the precession and BBH fits. The fit by BBH 
Model 1 is consistent with the fit by the precession, while BBH 
Model 2 cannot be reconciled with the precession fit at all. 
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Fig. 13. Two-dimensional path of C7. The precession fit and 
the fit by BBH Model 1 are presented. 
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Fig. 14. Two-dimensional path of C7 within 1 mas of the nu- 
cleus. The difference between the precession and BBH Model 
1 fit is due to the orbital motion of the black hole ejecting the 
jet. 

Table 5. Orbital period solutions for Va = 0.05 c-0.15 c 



Parameter 


Value 


V. [c] 


0.050 


0.075 


0.078 0.100 


0.125 


0.150 


To [yrs] 


770 


505 


480 375 


290 


250 



same range of Va as in Table|3 The fit to the optical light curve 
is shown in Fig.^] It does not depend on the actual value of 
To because ToV^' = const, and thus the same shape of the 
lightcurve can be reproduced by a number of combinations of 
Va and T^. 

5.5. Determination of tine Alfven speed 

In order to reconcile the two families of VA-dependent solu- 
tions (the precession solutions in Table 0] and the orbital pe- 
riod solutions in Table |5}, the value of the Alfven speed must 
be constrained. The periodic changes of the position angle of 
component ejection plotted in Fig. can be used for this pur- 
pose. The Alfven wave perturbations, propagating at Va <k Vc 




1990.5 1991.0 1991.5 1992.0 1992.5 1993.0 

Observed Date 

Fig. 15. Optical variability in 3C345 in 1990-93. Solid Unes 
show the fits by the BBH model. Individual peaks result from 
the orbital motion in the BBH system. 

cannot account for such a short period. The observed period 
^'short ~ 9.5 yr must therefore reflect perturbations of the mag- 
netic field caused by the orbital motion in the BBH system. In 
this case, T^ = Pshort^cO.Smas, where 5c,o.5mas is the Doppler fac- 
tor of the beam plasma measured at the same location (0.5 mas 
from the nucleus) at which the oscillations of the ejection angle 
are registered. The kinematic fit yields 5c,o.5mas = 50, which re- 
sults in the orbital period Pg = 480 yr. The corresponding value 
of the Alfven speed is then 

Va = 0.078 c . 

It should be noted that the Doppler factor 5c,o.5mas = 50 is sim- 
ilar to the Doppler factors implied by the observed statistics of 
superluminal motions near blazar cores (Jorstad et al. 12001 at . 

5.6. Parameters of tine binary system 

With the Alfven speed determined, the masses of the black 
holes can now be estimated from Tg and Tp. Since the orbital 
separation oc (Mi + M2)'^^, a family of solutions can be pro- 
duced that satisfy To = 480 yr, Tp = 2570 yr and reproduce 
the observed optical lightcurve of 3C 345 and radio properties 
of C7 (the shapes of the optical and radio lightcurves depend 
on the mass ratio 1^12 - M1/M2). The resulting solutions are 
presented in Table|6lfor five selected values of Mi . 

In order to identify an acceptable solution for the black hole 
masses, stability conditions of the accretion disk must be con- 
sidered. The rotational timescale of a precessing accretion disk 
in a binary system is given by a characteristic period (Laskar & 
Robutei n^93t 

The disk in a BBH system will be stable over timescales com- 
parable with the typical timescales of existence of a powerful 
radio source in AGN, Tagn ^10^ yr, if two conditions are satis- 
fied: 

1. Ml > M2, where Mi is the mass of the black hole sur- 
rounded by the accretion disk. 
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Table 6. Black hole mass solutions for To - 480 yr 



Parameter 


Value 


Ml [lO^'Mo] 


2.0 


6.0 


7.1 


20.0 


60.0 


M2 [10** Mq] 


4.0 


6.6 


7.1 


12.3 


23.1 


ro [pc] 


0.25 


0.32 


0.33 


0.44 


0.60 


Tdisk lyrs] 


179 


229 


239 


314 


430 



Notes: r„ - orbital separation between the black holes, T^^^y - char- 
acteristic rotational period of the accretion disk around the primary 
black hole. 
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Fig. 16. Range of acceptable solutions for the mass M\ of the 
primary black hole in 3C345. The acceptable solutions ex- 
ist within the range 7.1 x lO'^M© < Mi < 7.5 x lO^M©. 
Asterisks mark the solution for yL\ 2 - 1, for which M\ - M2 - 
1 A X 1O^M0. Shaded areas cover the values of M\, To and ii\ 2 
for which the accretion disk cannot remain stable over the typ- 
ical timescale of existence of a powerful radio source in AGN, 

Tagn < lO^r. 

2. Tdisk < 7^0/2 

These conditions restrict the mass of the primary black hole 
to a range 7.1 x lO'^M© < Mi < 7.5 x lO'^M©. The lower limit 
corresponds to yUi 2 - 1, and the upper limit corresponds to 
^'disk - ^"0/2 - 240 yr The acceptable range of Mi is very 
narrow, and any value within this range can be adopted. Since 
/ii_2 does not change significantly over the acceptable range of 
Ml (//i,2 = 1.02 for Ml = 7.5 x lO^M©), ;Ui,2 = 1 can be 
adopted for the binary black hole mass solution, yielding 

Ml = M2 = 7.1 X lO^M© 

listed in Tables|2land|6l 

6. Discussion 

The BBH model applied to the quasar 3C345 ties together, 
within a single framework, the variability of the optical flux 
and the evolution of the superluminal feature C7 resulted from 
a powerful flare that occurred in 1991. The model reconstructs 
the physical conditions in the plasma ejected into the jet during 
the flare, and the dynamic properties of the binary system of 
supermassive black holes. 



The binary black hole scenario was previously applied to 
3C 345 by Caproni & Abraham J2004I I who attributed varia- 
tions of the ejection angle (similar to those plotted in Fig. 
to precession of the accretion disk. This approach yielded an 
uncomfortably short precession period of 2.5 < Tp < 3.8 yr, 
very high black hole masses (Mi = 4-5 x lO^M©, M2 = 3- 
4x lO'^M©), and extremely small orbital separation = 0.017- 
0.024 pc, which is » 100 /?g of the primary black hole. In such 
a tight binary, the effective accretion radius of the secondary, 
ra = 0. 034-0. 050pc, is larger than the orbital separation, and 
it is very difficult to maintain the accretion disk intact. Indeed, 
it can be shown (Ivanov, Igumenshchev & Novikov 1998^ that 
the mass loss rate 1112 due to a shock that is formed after each 
passage of the secondary through the accretion disk may be 
very large, and for the BBH parameters of Caproni & Abraham 
f2004), it can reach m2 ~ Mi( VTOa^^^r', where Mi is the rate 
of accretion on the primary black hole and or* - a/ 10"^ is the 
disc viscosity. For weak accretion rates or small disk viscosity, 
1112 can even become larger than M\ , which would most likely 
lead to a subsequent destruction of the accretion disk. 

Another argument against associating the ejection angle 
variability with the disk precession is presented by the observed 
long-term kinematic changes in the jet. The apparent acceler- 
ations plotted in Fig.|6land the long-term trend present in the 
ejection angle changes plotted in Fig. indicate the presence 
of a periodic process on timescales of ~ 2000 yr. This observed 
periodicity, though very difficult to establish, agrees well with 
the precession period recovered by fitting the BBH model to 
the trajectory and flux evolution of C7. 

The trajectory of C7 (and indeed the trajectories of sev- 
eral other superluminal features monitored closely in the jet 
of 3C345; Lobanov & Zensus 1996 L96 LZ9j is clearly 
non-ballistic on timescales spanning almost 10 years in the 
observer's frame (and more than 100 years in the comoving 
frame). This is a clear observational argument against a very 
short precession period in 3C 345. The model developed in this 
paper attributes non-ballistic trajectories of the jet features to 
Alfven perturbations of the magnetic field, possibly modulated 
by the orbital motion in the system. 

The sum of the black hole masses (Mi_2 = 1.4 x 1O''M0) 
estimated in our model agrees well with the mass Mbh = (4.0 ± 
2.4) + IO^Mq obtained from measurements of nuclear opacity 
and magnetic field strength in 3C345 (Lobanov 1998^. Gu, 
Cao & Jiang ( 2001 1 report a higher value of Mbh = 7.6x lO^M© 
for 3C 345, from the properties of the broad line region. This 
estimate, however, depend strongly on the monochromatic lu- 
minosity (Mbh . ) for which the errors can be as large 

^ ^ ™ 5100A^ __, = 

as 30% even for the nearby objects (Kaspi et al. |2000). This 
error could be even larger for weaker objects (3C 345 is about 
100 times weaker in the optical than an average object from the 
sample of Kaspi et al. 2000), and overestimating L}^^^^^ by only 
50% would be bring the Mbh down to a value similar to the one 
derived in this paper 

Finally, let us consider the periodicity of the flares and com- 
ponent ejections in 3C345. A ~ 4yr period of the flaring ac- 
tivity is inferred from the observed variability of the radio nu- 
cleus (LZ99l, and it is related with the ejections of superlu- 
minal features (see Fig. |3}- Since this period is much shorter 
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Fig. 17. Characteristic timescales of disk activity in 3C 345, 
compared to the quasi-period of the nuclear flares and to the or- 
bital and characteristic disk rotation periods). The disk thermal 
instability operates at f ~ 4yr at distances of ~ 20 Rc- The me- 
chanical (sound wave) instability reaches the same timescale at 
~ 200 7?G, where it becomes comparable with the orbital and 
drift (accretion rate variations) timescales of the disk. A com- 
bination of these effects may cause the flare activity observed 
in3C345. 



than Tp - 2570 yr. To — 480 yr and Tdisk — 240 yr obtained in 
our model, there should be some other periodic process affect- 
ing the accretion disk on timescales of ~ 4yr. FigureFTTlshows 
four different characteristic disk timescales in 3C 345, for the 
black hole and disk parameters determined in Sect. |4] These 
timescales can be compared to the periodicity of the nuclear 
flares. It appears that several processes occurring in the disk 
at radial distances of 20-200 Rc may have similar characteris- 
tic timescales. Thermal instability at ~ 20 Rc and mechanical 
(sound wave) instability at ~ 200 Rc will develop at timescales 
of ~ 4 yr. In addition to this, the orbital and drift (accretion rate 
variation) timescales approach ~ 4yr at ~ 200 7?g- A combina- 
tion of these effects can cause the observed flare activity. 

Another possible explanation for the observed flares is an 
intermediate or stellar mass black hole or even a neutron star 
orbiting the primary black hole at ~ 800-1 200 Rc at a large 
angle with respect to the accretion disk plane. Close approaches 
to, or even passages through, the accretion disk should lead to 
formation of shock waves in the disk (Ivanov, Igumenshchev & 
Novikov 1998 ; Subr & Karas 1999 1, which in turn would cause 
an increase in the accretion rate resulting in a nuclear flare and 
injection of a dense plasma cloud into the jet. 

The periodicity of the flares may also be caused by a tur- 
bulence in the jet. Formation of a turbulent layer close to the 
beam-jet boundary can enhance the pair creation rate in the 
beam, which would be observed as the appearance, and propa- 
gation, of a bright feature embedded in the flow. 



7. Conclusions 

The model developed in this paper explains the variability of 
the optical flux and the kinematic and emission properties of 
the feature C7 observed in the relativistic jet of 3C 345 after 
a strong nuclear flare. The dynamic properties of the system 
are described in the framework of the orbital motion in a bi- 
nary system of supermassive black holes and the precession of 
an accretion disk around the primary black hole. The emission 
and physical properties of relativistic plasma produced during 
the flare are treated within the two-fluid approach postulating a 
highly relativistic e"^ beam propagating inside an e^p jet. The 
e^p jet transports a major fraction of the kinetic power associ- 
ated with the outflow from the nucleus of 3C 345, stabilizes the 

beam against plasma instabilities and shields it from inter- 
action with the external medium. 

Determination of the properties of the jet and the binary 
system in 3C345 is a difficult process involving two major 
steps. First, the precession characteristics of the accretion disk 
are determined by fitting the kinematic and flux density evolu- 
tion of C7. This yields a precession period Tj, = 2570 yr, which 
is in a good agreement with the periodicities inferred from the 
changes of kinematic properties of individual jet components 
and variations of the ejecting position angle in 3C 345. 

In the second step, the orbital and physical properties of 
the binary system are recovered by fitting the observed opti- 
cal variability. This yields an orbital period of Tp - 480. The 
consistency between the two fits is ensured by requiring the fit 
to the optical light curve to reproduce also the kinematic and 
flux density evolution of C7. The consistency argument con- 
strains the properties of radio and optically emitting plasma, 
with Fniin.opt = 12 and T^inj^d = 4. 

After ensuring consistency between the two fits, they are 
combined in order to constrain the Alfven speed Va = 0.078 c, 
and make it popsible to recover the physical properties of the 
system. In the last step, the conditions of disk stability are em- 
ployed to determine the individual masses of the binary com- 
panions. This procedure yields Mi = M2 = 7.1 x lO'^Mo. 
Finally, the model offers an explanation for the observed 4- 
year quasi periodicity of nuclear flares in 3C 345. This quasi- 
periodicity arises naturally from the characteristic timescales 
of the accretion disk surrounding the primary black hole. The 
overall success that the model has in explaining a variety of 
the observational properties in 3C 345 argues in favor of this 
explanation. 
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